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A systematic comparison of three approximate self-interaction corrections (SICs), Perdew-Zunger 
SIC, Lundin-Eriksson SIC and extended Fermi- Amaldi SIC, is performed for a model Hamiltonian 
whose exact many-body solution and exact local-density approximation (LDA) are known. For each 
of the three proposals we compare its implementation only for the potential, only for the energy (i.e., 
a post-LDA evaluation of the SIC energy), to none of them {i.e., a standard LDA calculation) and 
to both. Each of the resulting 10 permutations of methodologies is applied to 420 Hubbard chains 
differing in size, particle number and interaction strength. A statistical analysis of the resulting data 
set reveals trends and permits to assess the performance of each methodology. Overall, but not in 
each individual case, a post-LDA application of Perdew-Zunger SIC emerges as the recommended 
methodology. 
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I. INTRODUCTION 

Density-functional theory (DFT)ii2ii^ is advancing at 
a rapid pace, driven by the demand of ever more accu- 
rate electronic-structure calculations. To meet this de- 
mand, ever more sophisticated density functionals are 
being constructed. A simultaneous, and often contradic- 
tory, demand is the applicability to larger and more com- 
plex systems, met by computational and methodological 
advances in the implementation of density functionals. 
In this context, the comparison of different approximate 
density functionals and different modes of their imple- 
mentation is of key importance. The present paper pro- 
vides such a comparison for the case of self-interaction 
corrected functionals. 

In Sec. |TT] we recall three different self-interaction cor- 
rections. In Sec, mil we discuss four different implemen- 
tations, selected specifically to be able to compare the 
effect of self-interaction corrections on the potential and 
on the energy, in an unbiased way. In Sec. IIVI we report 
results from a statistical analysis of hundreds of calcu- 
lations with each of the possible permutations of correc- 
tion and implementation, for an exactly solvable model 
system whose exact many-body ground-state energy is 
known. Sec. |V] contains our conclusions. 



II. SELF-INTERACTION ERROR AND 
SELF-INTERACTION CORRECTIONS 

One electron does not interact with itself. Simple as it 
is, this fact is at the heart of much trouble in many-body 
and electronic-structure theory. Traditional density func- 
tionals, such as the local density, local spin-density and 
generalized-gradient approximations (GGAs), as well as 
many more recently developed functionals, such as hy- 
brids and meta GGAs, do have a spurious self-interaction 
error (SIE), and as a consequence predict a nonzero in- 



teraction energy even for a single electron. The exact ex- 
change energy of one electron exactly cancels its Hartree 
energy, so that the Hartree -I- exchange energy of meth- 
ods employing exact exchange is self-interaction free, but 
the correlation energy remains subject to self interac- 
tion even in methods employing exact exchange. Con- 
sequences of the SIE of approximate functionals are mul- 
tifarious, and include wrong asymptotics of approximate 
exchange-correlation (xc) potentials, energy gaps, ion- 
ization energies, electron affinities and transition-metal 
magnetic moments. 

Already in 1934, Fermi and Amaldi (FA)'' realized that 
Thomas-Fermi calculations, which use a local-density ap- 
proximation for the noninteracting kinetic energy and 
approximate the full interaction energy by the Hartree 
energy. 
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suffer from a large SIE, and proposed a first self- 
interaction correction (SIC), of the form 
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EH[n]=EH[n]--EH[n]. (2) 



Equation ^ is clearly a very crude approximation, vi- 
olating, among other things, the requirement of size- 
consistency, but it does eliminate the SIE of Eh [n] for a 
one-electron system (TV = 1). 

After the development of modern density-functional 
theory, the question of how to eliminate the SIE posed 
itself anew. One way to carry over to modern Kohn- 
Sham DFT the simplicity of the FA approach is what 
we call extended-Fermi- Amaldi SIC (EFASIC), where the 
FA factor is applied not just to the Hartree energy, but 
to the sum of Hartree, exchange and correlation energies. 
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where E^^P'^°'^[ni,ni] is an approximation, e.g. LDA or 
GGA, to the exact exchange-correlation functional 

A very successful proposal for a more sophisticated 
SIC, made in 1981 by Perdew and Zunger (PZ)jl is 
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where no-(r) = X^r^'^-fe'^Ci")' '^feo-(r) = \ipka{r)\'^, and 
(/jfecr(r) are the Kohn-Sham (KS) orbitals. 

In PZSIC, the SIE is subtracted on a orbital by 
orbital basis. Evidently, for a one-electron density 
^^/^■f'^[n(i)] = -Enin'^^^. The PZSIC functional 
depends on the orbitals fkcr, so that its minimiza- 
tion with respect to the density must employ the 
optimized-effective potential (OEP) method or one of its 
simplifications f^iii^ In practice, however, the minimiza- 
tion is usually performed with respect to the orbital den- 
sities, and 

VsM^^ir) = Ve.t{r) + VH[n]ir) + <PP;-[n.,7i^](r) 

-«ffK.](r)-<?P;-K^,0](r)(5) 
=: ve^tir) + VH[n]{r) + v^,%'^'[n^,n^]ir) (6) 

is adopted as corresponding single-body potential. This 
potential is constructed from the density of all electrons. 

The PZSIC approach, simplified by minimizing with 
respect to orbital densities instead of total densities, has 
been very succesful in removing the one-electron SIE in 
solids 1^ Results from proper minimization with respect 
to the densities, by means of the OEP, are reported in, 
e.g., Rcfs. Iiollll for atomic systems. PZSIC has become 
so popular that frequently the abbreviation SIC is used 
as synonymous with PZSIC, but it has been repeatedly 
noted that the PZ proposal is not the only possibility for 
removing the SIEJ^ii^iiiii^ 

An innovative proposal for an alternative self- 
interaction correction was put forward in 2001 by Lundin 
and Eriksson (LE).^^ The LE proposal attempts to con- 
struct an effective potential that, when acting on the or- 
bital of one electron, is constructed only from the density 
of the other electrons. In practice, this is done by intro- 
ducing an orbital-dependent effective potential, Vs^ka-i'<^)i 
determined by subtracting from the full density n(r) the 
partial density of the orbital the potential acts on^^^ 

vit^'^(r) = ve^tir) + VH[n- nfe,](r) 

+v^,r''^[n,-nka,n^]{r) (7) 

= Vext{r) + VH[n]{r) - VH[nka]{r) 

+<f°^[n,-nfe,,n^](r) (8) 

= : Vextir) + VH[n]ir) + v^^^^J^[n^, n^](r). (9) 

In this way, the approximate effective potential acting 



a one-particle system this approach also correctly zeroes 
the interaction contribution to the effective potential. 

The LE proposal for a corrected effective potential is 
accompanied by a similar expression for the exchange- 
correlation energyfiSiii^ 

occ 

i?i'f^'''K,^i] = -E E ^«K.] 

k <T=T,i 
occ p 

+ E E dr^nUr)e7c'"''"'[na-nka,n^]. (10) 

We note that the change from the LESIC to the PZSIC 
approach can be affected by substituting 
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on a given orbital. 



.LESIC 



(r), is constructed only from 



the density arising from the other orbitals. Clearly, for 



(12) 
in the xc energy density and potential of former. For 
the Hartree potential, the substitution vnin — rika] ^ 
VH[n] — VHirikcr] is an identity. 

The LE energy expression (jlOp is constructed in anal- 
ogy to Eq. ([9]) , but its functional derivative is not the LE 
xc potential in Eq. ([9]). Rather, both are separate con- 
structions. This gives rise to a certain ambiguity in how 
to implement the LE approach. Sec. lIIII we therefore dis- 
cuss various alternative implementations of Eqs. ([9]) and 

m- 

In the original LE worki^ it was argued that the LE 
approach should be superior to the PZ approach because 
it removes the SIE of the potential, which, according 
to Ref. 'IS, remains in the PZ approach. It has been 
objected^^ that the LE proposal cannot be right because 
it would even correct the hypothetical exact functional, 
for which the PZ approach correctly reduces to zero. 

While we agree with this objection as a matter of prin- 
ciple, we feel that in practice the key issue is not only 
what the correction does to the hypothetical exact func- 
tional, but also what it does to actually available approx- 
imate functionals. If such functionals were consistently 
improved by a correction, few workers would refrain from 
using this correction in practice, only because it overcor- 
rects the hypothetical exact functional. 

Interestingly, there are hints in the literature that ap- 
proximate xc functionals do indeed benefit more from 
LESIC than from PZSIC. Novak et al.^^ compare the per- 
formance of LSDA, LSDA+PZSIC and LSDA+LESIC 
in the calculation of hyperfine parameters, and find 
that LESIC significantly improves agreement with exper- 
iment, relative to PZSIC. (See also Ref. 19 for a success- 
ful application of LESIC to hyperfine parameters.) In 
a separate study, Friis et al.— compared density distri- 
butions predicted for the Mg crystal with experimental 
electron diffraction data, and observed that LESIC pro- 
duces better core and valence densities than PZSIC. 



The available data are too limited, however, to already 
conclude that LESIC is definitely superior, in practice, 
to PZSIC, in particular since, as a matter of principle, 
LESIC cannot be correct. It therefore becomes an impor- 
tant task to systematically investigate the performance 
of both approaches. The present paper is a first step to- 
wards this task, providing a systematic analysis of the 
performance of PZSIC, LESIC and EFASIC for a weU 
controlled model system whose exact solution is known. 



III. IMPLEMENTATION 

Formally, the xc potential corresponding to a given 
approximation to E^c [n^ , n^] is 



JT rpapprox 
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We have just seen, however, that by construction the LE 
xc potential is not the derivative of the LE xc energy, but 
a separate construction. Strictly speaking, the same is 
true for most common implementations of PZSIC, where 
the implemented potential is the orbital derivative, and 
not the density derivative of the energy. This ambiguity 
suggests that a complete test of each proposal involves 
four different permutations of methodologies. 

Our notation for these is APPV-APPE, where APPV 
represents the approximation made for the potential en- 
tering the KS equations, and APPE the approximation 
made for the energy functional. The simplest possibility 
is thus LDA-LDA, which is a standard LDA calculation 
using the L(S)DA potential and energyi^ The second 
scheme, denoted LDA-SIC, consists of a self-consistent 
LDA calculation, followed by a single evaluation of the 
SIE corrected energy functional on the LDA densities 
and orbitals. This strategy is also known as post-LDA 
implementation of SIC. By construction, it tests the per- 
formance of the SIC for the xc energy. The inverse pos- 
sibility, which wc label SIC-LDA, is a self-consistent SIC 
calculation, followed by a single evaluation of the uncor- 
rected LDA energy functional on the corrected densities. 
This procedure tests the performance of the employed 
xc correction for the potentials. Finally, SIC-SIC refers 
to a calculation using the SIC under study both in the 
potential and the energy. Starting from LDA-LDA, each 
of the other three possibilities can be implemented for 
each of the three SICs described in the previous section, 
resulting in a total of ten combinations of methodologies. 

In the self-consistent implementation of the PZSIC po- 
tential obtained from the orbital derivative of the PZ en- 
ergy and of the LESIC potential, which are required in 
schemes SIC-LDA and SIC-SIC, we face the additional 
difficulty that these potentials are not common multi- 
plicative potentials for all KS orbitals, but are orbital- 
specific, i.e. different for each orbital they act on. Such 
orbital-specific potentials, Ufc[{<Pj"M}](r), are not proper 
Kohn-Sham potentials, although they are still density 
functionals. 



A common multiplicative potential for all orbitals can 
be generated from the PZ and LE energy functionals by 
means of the OEP, but in the case of LESIC the potential 
resulting from the energy functional (jlOp would have no 
relation to the LE potential ([9]), which, as wc have seen 
above, is a separate construction. Hence, implementation 
via the OEP would implicitly test the proposed energy 
correction PU)) . not the proposed potential correction ([5]), 
required in the schemes SIC-LDA and SIC-SIC. 

A way to directly compare the orbital-specific PZ and 
LE potentials is to average the xc potentials over all or- 
bitals, weighted with the contribution each orbital makes 
to the full spin density, according tc^^ 

-PZSIC 
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The result are average multiplicative potentials arising 
from the orbital-specific LE or PZ potential. Of course, 
some information is lost in the averaging, but since our 
aim is to compare LE and PZ among each other, the er- 
rors of the averaging procedure are of secondary impor- 
tance, as long as all functionals are treated in the same 
way. 

We note that in the case of the PZ correction the same 
averaged potential is also obtained by applying the OEP 
algorithm to the PZ xc functional, which in the first step 
leads to 
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and then making the Slater approximation^^ 
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for the functional derivative. In this sense, our SIC-LDA 
and SIC-SIC implementations of PZSIC are approximate 
OEP calculations. As explained above, the same is not 
true for LESIC, where the OEP would produce the local 
xc potential coming from the LE xc energy (jlOp . which 
is not the same as the LE xc potential in Eq. ([9]) or its 
average (IT51) . In the present context, it is the averag- 
ing prescription which must be employed to test the LE 
proposal for the potential, not the OEP. 



Finally, the question poses itself by what standard we 
are to measure the performance of a SIC. Comparison to 
experiment is hard because the SIE itself is not an ex- 
perimental observable, and analysis of fine differences be- 
tween computational approaches relative to experiment 
is complicated by effects such as relativity, finite temper- 
ature, zero-point vibrations, sample purity, etc. that are 
of no direct relevance to the calculation, but automat- 
ically included in the experimental value. Even on the 
purely computational side, the issue is thorny, because 
the approximate SIC is applied to a functional, such as 
LDA, that is itself approximate, and it is not always easy 
to disentangle the error of the functional from that of its 
correction, or from that of approximations such as the 
sphericalization of atomic densities often used in PZSIC 
calculations. Moreover, calculations for real systems in- 
volve basis sets or numerical meshes, whose adequateness 
must be carefully tested, and which may introduce fur- 
ther residual uncertainties in the calculations. 

For all these reasons, we perform our analysis using 
as theoretical laboratory a very well controlled model 
Hamiltonian: the one-dimensional Hubbard model 
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FIG. 1: Root-mean-square error, in percent, of the dif- 
ferent implementations of the Extended Fermi- Amaldi self- 
interaction correction (EFASIC), relative to the exact many- 
body ground-state energy. Ns = 42 systems are considered 
at each U/t. 



IV. RESULTS 



This model, which may be considered a special case of the 
PPP model popular in quantum chemistry,^? describes 
electrons on a one-dimensional lattice of size L, inter- 
acting with on-site interaction C/, and hopping from one 
site i to the next with amplitude t. Occupation of each 
site is limited to two particles, necessarily of opposite 
spin. Like-spin particles do not interact in the most com- 
mon form of the model. The basic Hohenberg-Kohn and 
Kohn-Sham theorems of DFT hold for this model, too, 
once the density n(r) is replaced by the on-site occupa- 
tion number— i2^ 



n(r) 






Y^n,,^Y.^clc,-)- (20) 



In terms of this variable, local-density and spin-density 
approximations for Hubbard chains and rings have been 
constructed )^i21iSSi2^ and can be employed in the usual 
way. 

What makes this model attractive as a theoretical 
laboratory for DFT is that its exact many-body solu- 
tion is known in the thermodynamic limit (i — > cxd, 
rii = n = const) r^ and an exact numerical solution is 
possible for small systems {L ~ 14 on a workstation, 
for arbitrary density distributions Ui). As a consequence 
of the availability of the exact solution for the infinite 
homogeneous system, the exact LDA is known, too.'^^ 
Moreover, a complete basis consists of two spin-orbitals 
per site, so that we can always work at the basis-set 
limit. Another welcome feature, which we will exploit 
below, is that the model allows to continuously vary the 
interaction strength both in the exact and in the DFT 
calculations, which simulates the behavior of weakly and 
strongly interacting systems in ab initio calculations. 



In this section we compare the three self-interaction 
corrections EFASIC, LESIC and PZSIC, in the LDA- 
LDA, LDA-SIC, SIC-LDA and SIC-SIC schemes, using 
the orbital (Slater-type) average of Eqs. P^ and psp to 
deal with the orbital-specific PZ and LE potentials in the 
schemes SIC-LDA and SIC-SIC. All calculations are done 
for finite chains of different lengths L, particle numbers 
N and interaction strengths U. Variation of these pa- 
rameters for each combination of methodology produces 
thousands of data, which we analyze statistically by re- 
porting the percentage root-mean-square (rms) error. 



rms = 100 
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(21) 



of each approximation relative to numerically exact 
many-body energies for the same set of model parame- 
ters. Ns is the number of systems in the data set, labelled 
by j. 

Our results are summarized in Figs. [1] to [3] and in Ta- 
ble |T1 Each of the three figures compares the rms error 
defined in Eq. ([21]) as a function of interaction strength 
U, ranging from the noninteracting system [7 = 0, over 
weakly U ~ 2t, strongly U ^ 6t to very strongly, 
U ^ lot, interacting systems. For each value of U/t, 
our sample consists of Ng — 42 systems of different 
sizes, particle number, and densities (L = 3,4, 5. .14 and 
iV = 2,4,6, ..12, with7V/L< 1). 

We can analyze these figures from two different points 
of view. First, comparing the three corrections, we 
find that PZSIC consistently performs much better than 
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FIG. 2: Root-mean-square error, in percent, of the differ- 
ent implementations of the Lundin-Eriksson self-interaction 
correction (LESIC) , relative to the exact many-body ground- 
state energy. Ns — 42 systems are considered at each U/t. 



FIG. 3: Root-mean-square error, in percent, of the different 
implementations of the Perdew-Zunger self-interaction correc- 
tion (PZSIC), relative to the exact many-body ground-state 
energy. Ns — 42 systems are considered at each U/t. 



LESIC or EFASIC. The scale on the vertical axis of the 
three figures is not the same, to permit a better visualiza- 
tion of the differences between the various implementa- 
tion schemes, but the values of the rms error are directly 
comparable, and leave little doubt about the superiority 
of PZSIC for these systems. Somewhat surprisingly, the 
very simple EFASIC still performs a little better than the 
more sophisticated LESIC. Of course, our conclusions are 
based only on ground-state energies, but since the differ- 
ence in performance for the energy is quite large, and 
since many other observables are obtained from energies, 
we expect that PZSIC should be superior to the others 
in most applications. 

Second, comparing modes of implementation, we find 
that they fall in two groups, distinguished by the energy 
functional: The schemes LDA-LDA and SIC-LDA always 
give similar results, as do the schemes LDA-SIC and SIC- 
SIC. The two schemes in the first group evaluate the LDA 
energy functional on self-consistent (LDA or SIC) densi- 
ties, while the two schemes in the second group evaluate 
the SIC energy functional on self-consistent (LDA or SIC) 
densities. The formation of these two groups, occuring 
in the same manner for each of the three corrections, 
shows us that what matters for the final result is the en- 
ergy functional evaluated once on self consistent densities, 
not the potential employed during the iterations to self- 
consistency. This observation gives support to common 
post-LDA implementations of complex functionals. 

For EFASIC and LESIC, the implementations in group 
two perform much worse than those in group one. Within 
group one, SIC-LDA does shghtly worse than LDA-LDA, 
while within group two SIC-SIC is slightly worse than 
LDA-SIC. For all U/t, a simple selfconsistent LDA cal- 
culation is better than a calculation employing EFASIC 
or LESIC for the potential, the energy, or both. 



For PZSIC, a crossover between the error curves occurs 
around U ^ 5t. Below this value, the picture is the same, 
qualitatively, as for LESIC and EFASIC: it is better to 
do a simple LDA calculation than to apply PZSIC in any 
implementation. Above U ~ 5t, the situation changes, 
and the schemes in group two (LDA-PZSIC and PZSIC- 
PZSIC) perform clearly better than those in group onc.'^l 
Hence, only PZSIC is a real improvement on the LDA, 
and this improvement manifests itself only for sufficiently 
strongly interacting systems. 

As the data show, for weakly interacting systems 
(where the physics is dominated by the kinetic energy 
term, and the particles are delocalized) LDA is still bet- 
ter than PZSIC. Of course, PZSIC removes the self- 
interaction error of the LDA also in this regime, but 
PZSIC itself is not exact, and thus introduces other er- 
rors. For weak interactions, it does not pay to correct a 
small error by a correction introducing larger ones. For 
stronger interactions, the particles start to localize, the 
SIE becomes more important, and its removal by PZSIC 
improves the total energies considerably. The separation 
in these two regimes is a very gratifying result to ob- 
tain, as f/ ~ 5i is known to mark the region where the 
Hubbard model crosses over from weakly to strongly in- 
teracting systems, and the electrons become increasingly 
localized. The performance of PZSIC is clearly correlated 
with this crossover. In fact, in retrospect, one could have 
predicted that the electrons in the Hubbard model start 
to localize around C/ w 5i simply by comparing the rela- 
tive performance of LDA and LDA-I-PZSIC as a function 
of U/t. 

To make this division in two distinct regimes clearer, 
we report in Table [J the rms error for each regime, i.e., 
sum not only over systems with different size and particle 
number, but also over those with interaction strengths 



TABLE I; Root-mean-square error for all weakly and all 
strongly interacting systems (columns 1 and 2) and for all 
investigated systems (column 3), for all ten permutations of 
methodologies investigated here. Entries in boldface are the 
best in each column (interaction regime). 



Method 


U/t = 1..5 


U/t = 6.. 10 


U/t = 1..10 


Ns 


210 


210 


420 


LDA-LDA 


1.6008 


4.2957 


3.2416 


LDA-EFASIC 


4.3713 


7.3541 


6.0494 


EFASIC-EFASIC 


4.4266 


7.4423 


6.1230 


EFASIC-LDA 


1.6551 


4.3827 


3.3126 


LDA-PZSIC 


2.6724 


3.0596 


2.8726 


PZSIC-PZSIC 


2.7471 


3.2693 


3.0189 


PZSIC-LDA 


1.6567 


4.4749 


3.3741 


LDA-LESIC 


7.3476 


13.9210 


11.1306 


LESIC-LESIC 


7.4749 


14.3770 


11.4580 


LESIC-LDA 


1.7020 


4.5960 


3.4655 



in the indicated range. The lowest rms error in each 
regime is printed in boldface. Clearly, for weakly inter- 
acting systems, LDA-LDA does best, but PZSIC-LDA 
and EFASIC-LDA are only marginally inferior because 
the SIC chosen for the potential does not matter much, 
compared to the effect of the energy functional. For 
strongly interacting systems, on the other hand, LDA- 
PZSIC wins, closely followed by PZSIC-PZSIC. OveraU, 
i.e. considering the rms error over all considered systems, 
regardless of the value of U/t, LDA-PZSIC is the best 
combination of correction and implementation among all 
methodologies tested here. 

We finally remark that the peculiar aspect of the LE 
approach, where the SIC potential is neither the density 
derivative nor the orbital derivative of the SIC energy, 
but a separate construction, does not seem to cause any 
additional features in the results. As we found above, the 
effect of the potential of any SIC is much less important 
than that of its energy, so the inconsistency between them 
in the LE approach apparently does not matter much. 



V. CONCLUSIONS 

We compared selfconsistent LDA calculations, three 
self-interaction corrections applied to the LDA, and three 
ways to implement each correction, for a total of 420 sys- 
tems, for each of which we also obtained, numerically, 
the exact many-body solution. Three clear trends emerge 
from the resulting 4200 approximate and 420 exact data: 
(i) For all three investigated corrections, the main effect 
comes from the energy, not the potential. Post-LDA im- 
plementations (correcting only the energy) and selfcon- 
sistent implementations (correcting the energy and the 
potential) give very similar results, -2^ while implementing 
only the correction to the potential has virtually no effect 
on the resulting ground-state energies, (ii) For weakly in- 
teracting systems, the benefit of an approximate SIC is 
overcompensated by the intrinsic errors of the approx- 
imation. Uncorrected LDA here works best, (iii) For 
more strongly interacting systems, correcting the SIE is 
important, but of the three investigated corrections only 
PZSIC systematically improves on the LDA. Both the 
simple EFASIC and the more sophisticated LESIC are 
much worse. 

Of course, all our conclusions may have been biased 
by the choice of the Hubbard model as theoretical labo- 
ratory. The disadvantages of such possible bias, on the 
other hand, are largely compensated by the availability 
of an exact many-body solution as a benchmark, of the 
exact LDA,?^ and by the absence of any problems due 
to finite basis sets or due to experimental uncertainties. 
It remains to be seen if the superiority of PZSIC re- 
mains unchallenged by extension of these investigations 
to other corrections, implementations and classes of sys- 
tems. Similarly, the important issue of the many-electron 
gjj j33,34 j-gquij-es separate investigation. 

Acknowledgments This work was supported by 
FAPESP and CNPq. KC thanks J. Terra and D. 
Guenzberger for bringing the LE approach to his atten- 
tion. We thank V. L. Campo Jr. for his exact diagonal- 
ization code for the Hubbard model. 



Electronic address: 'capelle@ifsc.usp.br' 
W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 
R. M. Dreizler and E. K. U. Gross, Density Functional 
Theory: An Approach to the Quantum Many-Body Prob- 
lem (Springer- Ver lag, Berlin, 1990). 

R. G. Parr and W. Yang, Density-Functional Theory of 
Atoms and Molecules (Oxford University Press, Oxford, 
1989). 

E. Fermi and E. Amaldi, Acad. Ital. Rome 6, 117 (1934). 
See P. W. Ayers, R. C. Morrison and R. G. Parr, Mol. 
Phys. 103, 2061 (2005) for a modern discussion. 
J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 
J. B. Krieger, Y. Li, and G. J. lafrate, Phys. Rev. A 45, 
101 (1992). 



S. Kiimmel and J. P. Perdew, Phys. Rev. B 68, 035103 

(2003). 

W. Yang and Q. Qu, Phys. Rev. Lett. 89, 143002 (2002). 

A. Svane, L. Petit, Z. Szotek, and W. M. Temmerman, 
Phys. Rev. B 76, 115116 (2007). M. Liiders, A. Ernst, M. 
Dane, Z. Szotek, A. Svane, D. Kodderitzsch, W. Hergert, 

B. L. Gyorffy and W. M. Temmerman, Phys. Rev. B 71, 
205109 (2005). L. Petit, A. Svane, W. M. Temmerman, 
and Z. Szotek Phys. Rev. B 63, 165107 (2001). P. Strange, 
A. Svane, W. M. Temmerman, S. Szotek and H. Winter, 
Nature 399, 756 (1999). 

J. Chen, J. B. Krieger, Y. Li and G. J. lafrate, Phys. Rev. 

A 54, 3939 (1996). 

M. Norman and D. D. Koelling, Phys. Rev. B 30, 5530 



(1984). 
^^ P. Cortona, Phys. Rev. A 34, 769 (1986). 
" H. -J. Unger, Phys. Lett. A 284, 124 (2001). 
^^ S. H. Vosko e L. Wilk, J. Phys. B: At. Mol. Phys. 16, 3687 

(1983). 
^^ U. Lundin and O. Eriksson, Int. J. Quantum Chem. 81, 

247 (2001). 
" P. Novak, J. Kunes, W. E. Pickett, W. Ku and F. R. Wag- 
ner, Phys. Rev. B 67, 140403 (2003). 
" D.-K. Seo, J. Chem. Phys. 125, 154105 (2006). 
^* J. P. Perdew, A. Ruzsinszky, J. Tao, V. N. Staroverov, G. 

E. Scuseria and G. I. Csonka, J. Chem. Phys. 123, 062201 

(2005). 
^^ D. Guenzburger, D. E. Ellis and J. Terra, Theo. Chem. 

Ace. 113, 191 (2005). 
2° J. Friis, G. K. H. Maden, F. K. Larsen, B. Jiang, K. 

Marthinsen and R. Holmestad, J. Chem. Phys. 119, 11359 

(2003). 
^^ It is understood that whenever the system has unpaired 

spins, LDA is replaced by LSDA, and that in the present 

paper SIC always refers to SIC as applied to L(S)DA. 
^^ J. C. Slater, Phys. Rev. 81, 385 (1951). 
2^ R. Pariser and R. G. Parr, J. Chem. Phys. 21, 767 (1953). 

J. A. Pople, Trans. Farad. Soc. 49, 1375 (1953). 
■^* O. Gunnarsson and K. Schonhammer, Phys. Rev. Lett. 56, 

1968 (1986). 
^^ K. Schonhammer, O. Gunnarsson, and R. M. Noack, Phys. 

Rev. B 52, 2504 (1995). 
2"^ N. A. Lima, M. F. Silva, L. N. Oliveira, and K. Capelle, 



Phys. Rev. Lett. 90, 146402 (2003). 

'^'^ N. A. Lima, L. N. Oliveira and K. Capelle, Europhys. Lett. 
60, p. 601 (2002) 

2® K. Capelle, N. A. Lima, M. F. Silva and L. N. Oliveira 
in The fundamentals of electron density, density matrix 
and density functional theory in atoms, molecules and the 
solid state, N.I. Gidopoulos and S. Wilson eds., Kluwer 
Series Progress in Theoretical Chemistry and Physics, p. 
145 (2003) [also available as eprint cond-mat/0209245 . 

2^ G. Xianlong, M. Polini, M. P. Tosi, V. L. Campo, Jr., K. 
Capelle, and M. Rigol, Phys. Rev. B 73, 165120 (2006). 

'^° E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

^^ By exact LDA we mean that there is no need for analytical 
interpolations such as the PZ81, VWN or PW92 forms of 
the ab initio LDA, since the energy of the homogeneous 
reference system can be calculated numerically exactly for 
arbitrary densities. 

^^ We speculate that the small advantage the post-LDA 
scheme LDA-PZSIC has over selfconsistent PZSIC-PZSIC 
is lost once the common multiplicative PZSIC potential 
used in the latter scheme is implemented by means of the 
full OEP instead of its Slater approximation . ' 

^^ P. Mori-Sanchez, A. J. Cohen and W. Yang, J. Chem. 
Phys. 125, 201102 (2006). 

^■^ A. Ruzsinszky, J. P. Perdew, G.I. Csonka, O.A. Vydrov, 
and G.E. Scuseria, J. Chem. Phys. 125, 194112 (2006); 
ibid 126, 104102 (2007). 



